On scaling laws in turbulent magnetohydrodynamic Rayleigh-Benard convection 
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We invoke the concepts of magnetic boundary layer and magnetic Rayleigh number and use 
the magnetic energy dissipation rates in the bulk and the boundary layers to derive some scaling 
laws expressing how Nusselt number depends on magnetic Rayleigh number, Prandtl number and 
magnetic Prandtl number for the simple case of turbulent magnetohydrodynamic Rayleigh-Benard 
convection in the presence of uniform vertical magnetic field. 
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Turbulent fluid convection is an unsolved problem hav- 
ing very wide applications in the study of convective pro- 
cesses in atmospheres, oceans, metallurgy etc. If the fluid 
is conducting and is acted upon by magnetic field then a 
theory for this magnetohydrodynamic fluid's convection 
hopes to explain, in the long run, the convection pro- 
cesses in planetary core, stellar interior and; some other 
important astrophysical, industrial and geophysical situ- 
ations. Complex nature of these realistic situations have 
prompted the researchers first to try to find a theory 
for the rather simpler problem of Rayleigh-Benard (RB) 
convection [1], probably also because doing experiment on 
it is quite feasible. 

Briefly speaking, RB convection in MHD fluid is inves- 
tigated in a set up where a magnetofluid of density p, 
kinematic viscosity v, conductivity a, thermal diffusiv- 
ity k, magnetic permeability p and isobaric thermal ex- 
pansion coefficient a is confined between two horizon- 
tal plates (conducting or non-conducting). Theoretically, 
the plates are considered to have infinite extent; in real- 
ity, this can be realised by making the thickness d of 
magnetofluid between them very small compared to the 
lateral extent. The plates are maintained at a constant 
relative temperature difference A. The constant acceler- 
ation due to gravity g is acting in the downward direction 
and the entire set up is acted upon by constant uniform 
vertical magnetic field B. Experiments [2, 3] on MHD 
RB convection has been done using gallium or mercury 
confined between two copper plates. However, whether 
the results of ref.[3] applies to this paper is a question 
because this paper is applicable for high Nusselt number 
flows unlike the reference. 

If u(x, y, z, t), p(x, y, x, t), T(x, y, z, t) and b(x, y, z, t) be 
the velocity field, the kinematic pressure field (containing 
also the external divergence-free force terms), the tem- 
perature field and the magnetic field perturbation respec- 
tively, then the equations describing the convection under 
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the Boussinesq's approximation are: 

d t u t + UjdjUi = —dtp + vdjdjiii + agTS i3 
+ —Bjdjbi 



(1) 



dtbi + Ujdjbi = BjdjUi + bjdjUi H djdjbi (2) 

op 



d,T- 



,9,-T = ndidiT 



(3) 



The second order terms have been purposefully retained 
in these three equations. The boundary conditions at the 
two horizontal rigid plates, assumed perfectly conducting 
(electrically), are: a) u — at z — and z = d, b) T — 
±A/2 onz = and d respectively and, c) h z — Q on the 
plates. As the convection becomes turbulent, assuming 
the existence of wind of turbulence i.e., there exists a 
mean large scale velocity U that stirs the magnetofluid in 
the bulk, we can define a Reynolds number Re along with 
other important non-dimensional parameters. These are 
listed below: 

Reynolds number: Re = Ud/v 
Rayleigh number: Ra — gaAd 3 /nv 
Chandrasekhar number: Q — aB 2 d 2 / pv 
Prandtl number: Pr = i//k 
Magnetic Prandtl number: Pm = pov 
Nusselt number: Nu = (u z T - K,d 3 T) / (kA / d) 

where the overline denotes average over any horizontal 
plane. By the way, in this paper we shall use angular 
brackets to denote volume average. 

It has been conjectured[4] that there are various regimes 
in MHD RB convection wherein the scaling laws connect- 
ing Nu with Ra and Q are different; this fact has been 
verified later by experiments [2]. We shall take the re- 
sults of these experiments at their face values although 
some deeper meaning behind them must exist; most im- 
portantly we shall make the following two observations 
from the results of that paper: 

(i) The Nu mostly depends on some power of the ra- 
tio Ra/Q. This ratio basically is a sort of Rayleigh 
number — which we shall call magnetic Rayleigh 
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number Rb — constructed with magnetic viscos- 
ity aB 2 d 2 /p that essentially is the manifestation of 
Joule damping present in the magnctofluid. 

(ii) It is surprising that Nu — 1, rather than Nu, seems 
to be proportional to some powers of Rb. 

Of late, a paper [5] on turbulent RB convection in or- 
dinary fluids has presented a unifying theory on how 
Nu depends on Ra by constructing four regimes each of 
which is defined depending on what (the boundary layers 
or the bulk of the fluid) the dominant contributors to the 
kinetic and the thermal dissipation rates are. We shall 
closely follow this very paper extending the arguments 
and the assumptions to the turbulent RB convection in 
magnctofluid in the presence of uniform vertical magnetic 
field keeping in mind the observations (i) and (ii) made 
above. 

To begin with, we define e u = u((diUj)(diUj)) , £t = 
K((diT){diT)) and e h = (^ap)- 1 ((d^)^)) . Using 
equations (1) and (2), and the definition of the Nu, we 
arrive at: 



v s Ra 



(Nu - 1) 



d 4 p r 2 

Similarly, usage of equation (3) yields: 
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K 



Nu 



(4) 



(5) 



We shall now clearly state the main assumptions in- 
volved. The very first and most important assumption 
of this paper is: 



v z Ra 
v 3 Ra 



(Nu - 1) 
(Nu - 1) 



(6) 
(7) 



where (and henceforth, '~' means 'scales as'). This ex- 
plains (or is explained by) the observation (ii) above. 
Secondly, we assume that there is a magnetic boundary 
layer (BL) of thickness 5 b where the perturbation mag- 
netic field is strongly affected by the magnetic diffusivity 
so that it grows to some definite value at the top of the 
layer from zero on the surface of the plate. Thirdly, to 
go with the observation (i) above we shall churn this the- 
ory in such a manner that Nu depends on Rb and not 
on Ra explicitly. Fourthly, we shall confine ourselves to 
the laminar boundary layers only though the bulk of the 
magnetofluid is in turbulent state. So the thickness of 
the Blasius type kinetic BL 6 U and the thermal bound- 
ary layer St are given by 6 U ~ d/^J~Re and 8t ~ dNu^ 1 
respectively. Lastly, we shall confine our attention to 
the relevant cases where convection dominates conduc- 
tion i.e., Nu 3> 1. 

Armed with these hypotheses, we propose that there may 
be following eight possible different regimes depending 
on whether the bulk or the respective boundary layers is 
chief contributor to dissipations: 



(A) : e h (bl), e T (bl), e u (bl). 

(B) : e b (bl), e T (bl), e u (bulk). 

(C) : £ 6 (60, e T (bulk), e u (bl). 

(D) : Eb(bl), er(bulk), e u (bulk). 

(E) : e b (bulk), e T (bl), e u (bl). 

(F) : Eb(bulk), £r(bl), £ u (bulk). 

(G) : Eb(bulk), £x(bulk), E u (bl). 

(H) : £ b (bulk), e T (bulk), E u (bulk). 

Here, '(£>/)' and '(fru/fc)' respectively says that the respec- 
tive boundary layer and the bulk is dominant contributor 
to the dissipation concerned. These regimes can be pre- 
pared by proper choice of Ra, Pr and Pm. For example, 
for large Ra regime (H) is possible; for small Pr and large 
Pm regime (F) is expected and so on. 
Before proceeding further, let us estimate the width 
of the magnetic boundary layer. For Pm <C 1 (i.e., 
S u <C 6b), we may neglect the first and the second terms 
in the R.H.S. of the equation (2) and under steady con- 
vection one is left to compare the second term of L.H.S. 
to the only unneglected term in the R.H.S: : 



Ujdjhl r 
B II 



1 

ua 



U 



uo$l 
d 



RePm 



(8) 
(9) 
(10) 



Note that while estimating gradients we are using the 
magnitude of the uniform external magnetic field B 
rather than the typical magnetic field fluctuations. This 
also is an assumption in analogy with the standard as- 
sumption made in ordinary fluid RB convection where 
while estimating the temperature gradients A is used 
freely. For Pm 3> 1 (i.e., S u 3> 6b), we must balance 
the first and the third term in the R.H.S. of the equation 
(2): 



\BjdjUi\ 
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ua 



djdjbi 
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6 U 
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po6l 
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ReiPm^ 



(11) 
(12) 

(13) 



where the gradient of the velocity field in the kinematic 
boundary layer has been approximated by a linearly in- 
creasing profile. 

Now under the assumptions made, following can be easily 
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derived 



For Pr < 1, e T {bulk) 



e u (bulk) ~ ^- = ^Pe 3 (14) 



For Pr > 1, e T {bulk) 



UA 2 S T 
d S„ 



For Pm < 1, e b {bulk) 



'AV PrPef 

P 2 17 z/ 3 QPe 
2^7 ~ rf 1 Pm 



t, 1 ^ ,m P5 T ^ 3 QPef 

ForP m »l, ^fc).—--.- — 



(16) 
(17) 
(18) 
(19) 



Though we are interested in the limit of large and small 
Prandtl numbers (magnetic and non-magnetic), they can 
not be too large or too small, for, then cither the con- 
vective flow will be suppressed or the thermal and the 
magnetic diffusivities will be so high that the respective 
effects will be the dominance of thermal conduction and 
too quick decay of magnetic fields to have any magnetic 
effect on the fluid. In writing the equation (17), it has 
been kept in mind that the magnetic energy per unit mass 
B 2 /2pp is cascaded down the spatial scales by the turbu- 
lent eddies until it is dissipated at the lowest scale (the 
turbulent condition is taken to be stationary with contin- 
uous feeding of magnetic energy into the magnetofluid) 
and in the equation (19), the factor S u /d comes in to as- 
sert that only the kinematic boundary layers' volume is 
involved. Attempt to get similar relations for ex(bl) takes 
one to the expression (5). Therefore, alternate relations 
must be sought. Using the equation (3) and the argu- 
ments given in the Landau's text on fluid mechanics [6], 
one readily arrives at: 



ForPr<l, Nu~ Re^Pr* 
For Pr > 1, Nu^ReiPri 



(20) 
(21) 



Again, coming on the scaling expression for £&(W) we 
work out as follows: In the magnetic boundary layer if j 
is the current density, then from the relation V x b = p,j 
one can estimate j as B/Sbp,- This means that mean rate 
of dissipation of magnetic energy per unit mass j 2 /ap is 
given as: 



e b (bl) 



B 2 



upp? 5? 



(22) 



This equation alongwith with relations (10) and (13) 
gives £b(bl) for Pm <C 1 and Pm >• 1 respectively. 
Now we are prepared to derive a range of scaling laws. 
We shall demonstrate the line of attack for doing so by 
finding out a illustrative scaling law for the regime (A); 



TABLE I: For Pm«l and Pr < 1. (UD=Undetermined) 



Regime 


Exponent of Rb 


Exponent of Pr 


Exponent of Pm 


A 


i 

1 








B 










C 


UD 


UD 


UD 


D 


UD 


UD 


UD 


E 


1 


-1 


1 


F 


1 


-1 


1 


G 


UD 


UD 


UD 


H 


UD 


UD 


UD 



TABLE II: For Pm < 1 and Pr > 1. (UD=Undetermined) 



Regime 


Exponent of Rb 


Exponent of Pr 


Exponent of Pm 


A 


i 

1 


•2 

9 





B 


§ ' 







C 




' 





D 


I ' 


_2 





E 


1 


% 

1 


1 


F 


1 




1 


G 


3 


-4 


3 


H 


3 


-4 


3 



all other seven regimes can be treated similarly. For the 
regime (A) the dominant contributor to the magnetic en- 
ergy dissipation are the magnetic BLs. Deciding to work 
for Pm 3> 1 first and as Nu >• 1, we equate the corre- 
sponding relation (22) (using the expression (13) for 5b) 
to the expression (7) to get: 



-^RaNuPr- 2 



B 2 



Nu • 



1 

Rb 
1 Pr 2 
Rb~Tm~ 



Pr \ 
Pm) 



app 2 5 2 

2 / , \ 2 



Pe5 



(23) 
(24) 
(25) 



Again in the regime (A), the dominant contributors to 
the kinetic energy dissipation are the kinetic BLs so we 
equate relation (19) to relation (6) to get: 



Re? = RaNuPr" 



(26) 



Also, they are the BLs — thermal BLs — that are con- 
tributing to the thermal energy dissipation, so relation 
(20) is of importance (we are considering Pr < 1 to 
begin with) for this case. For other cases, where the con- 
tribution to the thermal energy dissipation comes from 
the bulk, the corresponding er(bulk) will have to be com- 
pared with equation (5) to obtain the desired relations. 
Substituting relation (20) into relation (26), we obtain: 

Re^Ra^Pr~i (27) 

which in turn when put into relation (26), yields: 

Nu~ RaiPr-i (28) 
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TABLE III: For Pm»l and Pr < 1. (UD=Undetermined) 



Regime 


Exponent of Rb 


Exponent of Pr 


Exponent of Pm 


A 


i 

2 


l 

f 


i 

2 


B 


1 

2 




1 

2 


C 


UD 


UD 


UD 


D 


UD 


UD 


UD 


E 


2 


5 

2 


3 


F 


2 


5 

2 


3 


G 


UD 


UD 


UD 


H 


UD 


UD 


UD 


LABLE IV: For Pm > 1 and Pr > 1. (UD=Undetermined) 


Regime 


Exponent of Rb 


Exponent of Pr 


Exponent of Pm 


A 


i 

1 - 


i 

? 


i 

? 


B 


2 


2 


2 


C 


1 


-1 


1 


D 


1 


-1 


1 


E 


2 


-3 


3 


F 


2 


-3 


3 


G 


UD 


UD 


UD 


H 


UD 


UD 


UD 



Now, in accordance with the observation (i), we eliminate 
Ra from relations (27) and (28) to write Re in terms of 
Nu and Pr. This is going to be the general strategy 
throughout. Thus, for this example we have: 

Re ~ Nu 2 Pr~ l (29) 

Putting relation (29) in relation (25), we arrive at the 
desired scaling law: 

Nu~ RbiPr-*Pm? (30) 



The benefit of this strategy of finding scaling laws is 
that the exponents of Pr and Pm are also predicted. 
Moreover, the method of classification of the regimes 
in the present-day experiments [2] differs from what 
has been done theoretically here. Hence, each of those 
regimes may consist of an overlapping region of some of 
the eight regimes proposed in this paper. This would 
suggest that a fit of the form Nu = cRa a /Q b (a, b, c are 
just numerical constants) usually done to represent the 
experimental results can in principle be an equally good 

fit like: Nu = Cl Rb dl + c 2 Rb d2 + c 3 Rb d * H ; this of 

course is not mathematically impossible. 
In closing, we list rest of the all mathematically possible 
scaling results (assuming Nu <~ Rb a Pr b Pm c ) in four 
tables — table-I to table-IV — without showing explicit 
derivation whose strategy, anyway, has already been 
clearly outlined; and hope that experiments and simula- 
tions would be done in near future to test the conjecture 
proposed in this paper and the exponents of Prandtl 
number and magnetic Prandtl number will be concen- 
trated upon more seriously. Only then one can say if 
the idea presented herein, which seems to have a flavour 
of being a mere translation of the Grossmann-Lohsc 
theory [5] of thermal flow in Rayleigh-Bcnard geometry 
for similar convection processes in magnetohydrody- 
namic flows, is valid or not. Wc remind the readers 
that it has not been possible for us to investigate the 
scaling laws when the laminar BLs turn turbulent: This, 
of course, should draw attention of the theoretical fluid 
dynamists. 
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